Collective excitations and marginal stability of quantum Ising spin glasses 
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We solve the Sherrington-Kirkpatrick (SK) model in a transverse field V deep in its quantum glass phase at 
zero temperature. We show that the glass phase is critical everywhere, exhibiting collective excitations with 
a gapless Ohmic spectral function. Using an effective potential approach, we interpret the latter as being due 
to disordered spin waves which behave as weakly coupled, underdamped harmonic oscillators. In the limit of 
small transverse field V the low frequency tail of the spectrum tends to a universal limit independent of T. 



Spin glasses are canonical representatives of a wide class 
of complex disordered systems where competing (e.g. ferro- 
versus antiferromagnetic) interactions induce frustration, sup- 
pressing the emergence of simple ordering patterns. Never- 
theless the interactions induce a phase transition from a disor- 
dered paramagnetic to a glassy "ordered" state at low temper- 
atures. The emerging glass phase features several remarkable 
properties, which have been studied in great detail in classical 
glasses: despite of long range correlations there is no regular 
spatial order. Furthermore, the free energy landscape is very 
rough containing a large number of local minima, separated 
by high barriers [ ]. The latter entails the breaking of ergod- 
icity and intriguing long time out-of-equilibrium phenomena. 
Classical glass states feature criticality [ , ] in the sense that 
spin-spin correlations are power law correlated, despite the 
absence of a broken continuous symmetry. In systems with 
long range interactions (e.g., in Coulomb glasses, or in the 
SK model) this criticality is reflected by a pseudogap in the 
distribution of local fields at low temperatures [- -, ]. 

It is important to understand how such frustrated, critical 
systems behave in the presence of quantum fluctuations, and 
how the latter influence the dynamics and the relevant low en- 
ergy excitations, both close to quantum glass transitions, as 
well as deep in quantum glass phases. Experimentally, these 
questions are directly relevant for quantum spin glasses such 
as the compound LiHoa;Yi_a;F4, where quantum fluctuations 
can be tuned by an external magnetic field | |. Coulomb- 
frustrated semiconductors close to a metal insulator transi- 
tion [-f, 1 0- ], flux-frustrated Josephson junction aiTays [ 1 4- 
18], proton glasses (systems with frustrated ferroelectric inter- 
actions) [ ' - ], many-body cavity QED [ '], etc. Strong 
enough quantum fluctuations induce tunneling below the bar- 
riers and restore ergodicity, producing a quantum glass tran- 
sition [26, 27], which has been the subject of many investiga- 
tions [20, 21, 28-31]. 

In the past many theoretical studies on quantum glasses 
have focused on short range interacting systems and the en- 
hanced relevance of Griffith effects on the low frequency dy- 
namics in short range spin glasses [ ]. However, several of 
the above mentioned condensed matter realizations of frus- 
trated systems feature longer ranged interactions and may be 
approached theoretically from the limit of mean field mod- 
els, such as quantum rotors [33], SU{N) Heisenberg spin 
glass [ ] and the transverse field Ising spin glass [20, 27, 35- 



]. For the mean field version of the latter, the behavior at 
the quantum phase transition is well understood [ ], however 
little is known about the quantum glass phase except in the 
vicinity of criticality. Based on a Landau expansion, the glass 
phase has been conjectured to be critical [ ] with a spectral 
function behaving as |a;| at low frequencies. Exact diagonal- 
ization in small systems p^] has indeed shown such a trend. 
The quantum glass phase is expected to have similar complex- 
ity as the classical spin glass with multi-valley free energy 
landscape [35, 36, 41], but a formal description has remained 
elusive. So far, tractable models could be solved exactly only 
in the limit of a large number of vector/rotor components, in 
which the complexity and criticality captured by full replica 
symmetry breaking is lost. The latter phenomena are however 
among the most intriguing aspects of quantum glasses with 
critical phases. In contrast, quantum systems with closer sim- 
ilarity to structural glasses often exhibit discontinuous quan- 
tum glass transitions and non-critical glass phases [42]. 




FIG. 1. Phase diagram of the mean field quantum spin glass. The 
deep quantum glass (red dot) exhibits gapless collective excitations 
with a low frequency spectrum independent of the transverse field F. 



The present paper provides the missing link between the 
quantum phase transition and the deep quantum glass phase. 
We focus on Ising spin glasses, which combine a continuous 
glass transition with a critical glass phase and non-trivial er- 
godicity breaking. This model acts as a prototype for many 
glasses with long range interacting, discrete degrees of free- 
dom, such as localized electrons, electric dipoles etc. More- 
over, it has been recently been pointed out that the mean field 
model may be realized faithfully in random laser cavities. 



2 



where a multitude of modes provide random long range cou- 
plings between trapped atoms, as in the Dicke model [ ]. We 
study the quantum glass phase of the Sherrington-Kirkpatrick 
model in a transverse field [ ], i.e., the mean field version of 
the Ising quantum spin glass, 

i<i i 

where and cr^ are Pauli operators. Every spin interacts 
with all the others, and the couplings Jij are random Gaus- 
sian variables of zero mean and variance J"^ /N . The quantum 



The saddle point values of the off-diagonal Qab are time in- 
dependent [4,")] and have an ultrametric structure parameter- 
ized by the monotonous function [36, 41, 46] < q{x) < 
1, with X £ [0, 1] measuring the distance between replica 
in phase space, and thus being a proxy of a time scale 
in the aging regime. Rr-r' is the connected part of the 
replica-diagonal spin correlator Qaa{T, t'), which tends to the 
Edwards-Anderson parameter ^ea = 9(1) at large time sepa- 
rations at T = 0. Assuming a continuous function q{x), the 
self-consistency problem is equivalent to the solution of the 
diffusion-like equations [47, 48]: 

q(x\ ff f 
m[y,x) = — [m"{y,x) + 2Pxm{y,x) ■w!{y,x)] , 

P{y. x) = -^ [P"iy. ^) - 2/32; im{y, x) P{y, x))'] , 

q{x) = j dy P{y,x)m^{y,x), (3) 

where dots and primes denote derivatives with respect to x and 
y, respectively. P{y, x) is the distribution of frozen exchange 
fields y, averaged over time scales corresponding to phase 
space distances x, with P{y,0) — 6{y). Likewise, m{y,x) 
is the magnetization of a spin in the presence of a frozen field 
y on that time scale. Short-time observables are described 
by X = 1, of which the local field distribution P{y, 1) and 
the spin-spin correlator R^- will be of particular interest. The 
difference to the classical problem lies in the modified set of 
boundary conditions which read: m{y,x = 1) = {<J^)s{y) 



fluctuations are tuned by the transverse field T. The phase dia- 
gram and the deep quantum glass regime of particular interest 
to us are shown in Fig. 1 . In the classical limit F = a glass 
transition takes place at = J- As has long been known, it is 
connected by a line of continuous transitions to the quantum 
glass ti-ansition at Fc(T = 0) w 1.52 J [39, 43]. Below we 
use units with J = I, and restore J occasionally for clarity. 

Replica solution - We first solve the model by the replica 
approach, and then interpret its features with the physically 
more transparent effective potential (TAP) method [42, 44]. 
The disorder average of the free energy is carried out using 
the replica trick following Ref. [ ], reducing the problem to 
an effective, self-consistent single spin model: 



(2) 



I 

with the local action 

Siy) = ^ j j^Tdr' alRr-r'K' + driya'^+Va:). 

(4) 

The system of equations is closed by the self-consistency re- 
lation 

Rr-r' = I P(y,l)(T<a^,)5(,) -9(1). (5) 

Solution atT — Q- We now focus on the T = behavior 
Like at any quantum critical point, the gap closes at the tran- 
sition [ ]. However, it was found within replica symmetric 
Landau theory [ ] that the gap remains closed in its vicinity, 
with an Ohmic spectrum at small frequencies. Similar behav- 
ior was found in the analysis of a rotor model, where in the 
limit of M — >■ c» components [ ] the replica symmetry is 
not broken. Here we show that in the Ising case (believed to 
be the M = 1 limit of the rotor model) full replica symme- 
try occurs, and that the latter guarantees that the gaplessness 
extends into the whole glass phase. We point out that this phe- 
nomenology contrasts with that of the exactly solvable model 
SU(A^ — ?> oo) Heisenberg spin glass [ , ], which exhibits 
a random first order transition with distinct dynamic freezing 
and thermodynamic glass transitions. Its thermodynamically 
dominant states are gaped and are thus very different from the 
states obtained in the Ising limit — > 1 analyzed here. 

The spectral function is encoded in the Fourier transform 
Ruj of the average spin correlator Rr, which we analyze fol- 
lowing Miller and Huse [ ]. Representing the action S{y) 



^ (M! ^ + 1 /r drdr'Ql^ir, r') - lim 1 log Tr Te^^- 



Ml 



5eff= ffdrdr' 5]ga.<(r)a,^(r') + ^5]Qaa(r,r')<(r)<(r') + Vj^ f dr a^ir), 

Qaa{r,T') - (Ta,-(r)<(T')) = Cr-r' = Rr-r'+q{l). 



La<b 

Qab - («) q{x)- 
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with fermions, we can expand the spin-spin correlator into a 
power series in 



(6) 



where n;^(y) is the proper polarizabihty [4')], itself a func- 
tional of Ruj. However, 11 remains analytic at small uj, IIi^ ^ 
Ho—auj^, even when turns non-analytic in the glass phase. 
The latter is a consequence of the marginal stability, 



dyP{y,l){m')^{y,l) = 1, 



(7) 



which is implied by Eq. (3). Indeed, with the small frequency 
expansion R^ ^ Rq + SR^^, and noting that m'(y, 1) — 
Xu^oiy) as well as Xuiv) = Xoiv) + Xo(y)^^c. + 0(^2), 



Eqs. (5-7) require that 



This implies the non- 



analyticity Ru^o = Ro — B\uj\ in the low frequency cor- 
relator. Upon analytic continuation, the spin spectral function 



Buj 



ttN 



-iS 



is found to be always gapless and Ohmic, even deep in the 
Ising glass phase. Remarkably, as we will derive below, B 
becomes F-independent as F — > 0. The marginal stability and 
the related gaplessness of the replicon mode are both natu- 
ral by-products of continuous replica symmetry breaking in 
mean field glasses. They replace the role of Goldstone modes 
in systems with a broken continuous symmetry. Physically, 
the gapless excitations (8) arise due to the presence of many 
metastable states which are close to the ground state and con- 
nected to it via nearly flat directions in the energy landscape. 
The transverse field hybridizes the nearby classical states, 
which gives rise to very low-lying collective spin excitations. 

Deep glass phase - To obtain quantitative results we have 
solved the full self-consistency problem, focusing on the deep 
quantum glass phase where J ^ F ^ T 0. The limit 
T — > is taken by replacing the variable x € [0, 1] by 
/3cfr = /3x ^ [0, oo], which has the interpretation of an inverse 
effective temperature in the aging dynamics [50]. In the limit 
F <C J, the flow of Eqs. (3) is attracted to a scaling regime 
where dq/d(3^g = c(/3off)//?off' M^^v) Ml^osy), and 
P{x, y) /3cffp(/3offy), which holds for 1/ J < /3off < 1/F 
and y J. Here, c(/3cff) — >■ 0.411, and rh and p are the 
same fixed point functions that appear in the classical low T 
limit [:i 1]. In particular. 



p(/3cff2/) 



a|/3cff2;| 
const. 



1 < Wcsyl 
\Pcsy\ < 1, 



(9) 



which displays a linear pseudogap with slope a — 0.301 in the 
distribution of frozen fields, smeared on the scale y ^ l//3cff- 
For /3cff ^ 1/F, the overlap ceases to scale. The derivative 
drops [c(/3off) 0, cf. Fig. 2], g(/?eff) reaches the constant 
value qea, and the flow of P and m freeze to their short time 
forms. This yields the physical result that the frozen field 




y6eff/(/?eff+l/r) 

FIG. 2. Rescaled derivative of the overlap function, c = 
l3^gdq/dPcS as a function of it = /3cff /(/?cft' + 1 /F), extracted from 
the full solution for finite F (solid curves). The dashed curve is the 
solution of the asymptotically exact scaling anzats for F/Fc — > 
[ ]. The transverse field F takes values (bottom to top) F/Fc = 
0.05, 0.08, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. Convergence to the solution is 
slower [5 1,52] and numerical instabilities are more pronounced for 
lower values of F, F < O.lFc producing rougher curves. The scaling 
regime obtains at /3cff ^ 1/F (c — )■ 0.411). It reflects the ultramet- 
ric structure of phase space and an ensuing self-similar structure of 
dynamics deep in the aging regime [13, 51]. q{x) reaches its plateau 
value (jEA at = Xc/T ^ 0.5/F where c — >■ 0. 




FIG. 3. Opening of the linear pseudogap in the distribution of frozen 
fields in the glassy ground state at T = 0, P{y, 1), deeper and deeper 
into the quantum glass (F/Fc — >■ from top to bottom, same values 
as in Fig. 2). The inset shows the rescaled distribution P{y — zT) /F 
(solid lines), compared to the asymptotically exact scaling function 
p{z) (dashed line) [52]. 



distribution in typical glass states, P{y, 1) has a pseudogap, 
which is smeared on the scale F, as one may expect from sta- 
bility arguments. The evolution of the pseudogap upon enter- 
ing deeper into the quantum glass is shown in Fig. 3. 

In the limit F ^ J, F is the only relevant dynamic energy 
scale, while J merely determines the width of the distribution 
of frozen fields. We demonstrated this by analyzing the natu- 
ral scaling anzats: (i) for /3cff ^ 1/J, c(/3off) = c{u) is only 
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a function of u = /3off/(/3off + l/F); (ii) R{uj) = rf{uj/T) 
and (Hi) P{y,l) = Tpiy/T) [and P = {T/u)p{yu/T,u), 
m — rh{yu/r, u)]. One then verifies that c. f,p and m satisfy 
self-consistent equations which are independent of F <C r^. 
They admit a solution to which the full solutions of Eqs. (3)- 
(5) converge as F — > 0, as shown in Figs. 2 and 3. This so- 
lution implies a non-trivial scaling of the dynamic properties 
of the glass; in particular it entails the remarkable result that 
the coefficient B of the Ohmic spectral function (8) tends to a 
constant B w 0.59/ as F ^ 0. 

Physical interpretation - We now interpret the above re- 
sults physically, based on the effective potential (TAP) ap- 
proach [44]. The F-independence of the low frequency tail 
of the spectral function, predicted by the replica solution for 
the deep glass phase, will serve us as a benchmark. Using 
the method of Refs. [42, 53] we construct the Gibbs potential 
G[m, Cr] describing the free energy of the system constrained 
to have a magnetization pattern {mi] and global autocorrela- 
tion function Cr = l/^E,('7f (T)c^f (0)) (at T = 0): 



\NP J dT{Cr-q^Af 



(10) 



Here ^ea = l/^^Ei "^f is the free energy of single, 

constrained spin. The magnetizations of local minima nrii ~ 
(cr?) are computed self-consistently via SG/Snii = 0: 



dCo 



drrij 



Jijirij + mi x° 



0, 



(11) 



where x" = dT{Cr - Qea) = l/^EjXi is the static 
susceptibility. However, for quantum problems Eq. (11) is 
not closed since Go depends on the global autocorrelation 
function G(t), which has to be evaluated self-consistently, 
too [ ]. Since this exact formalism is too involved to yield 
direct physical insight, we approximate the static susceptibil- 
ities Xi and the local functional Gq by those of single spins, 
whose magnetization is constrained by an auxiliary static 
field: 



Go[m,] = -F(l 



2N1/2 



Xi = 



(1 - m|)3/2 



(12) 



This approximation is similar but not identical to the "static 
approximation" employed in replica approaches to quantum 
spin glasses [ ]. It overestimates the susceptibility to longi- 
tudinal fields, enhancing the stability of the glass. However, 
it reproduces qualitatively the results of the rigorous replica 
theory, furnishing a complementary physical picture. 

Collective excitations in a local minimum are governed by 
the curvature of the energy landscape, i.e. by the Hessian 



5'G 
SmiSmj 



1 

IXi 



J'ix) 



where (x*^) = ^/N^-Xi- The replica theory assures that 
the glass phase is marginal. Here this translates into a gapless 



spectrum of eigenvalues of Hij, which requires [54] 

J'{x')^i- 



(13) 



This is the natural analog of Eq. (7). Under this condition, the 
density of eigenvalues of the Hessian A is given by 



P(A) 



AF 



J2 



for A, - F < J. 



To establish the link with the spin spectral density (8) we 
interpret the low energy normal modes of H as weakly inter- 
acting [ ] harmonic oscillators with spring constant A, an ef- 
fective mass scaling as M ^ 1/F, and thus an eigenfrequency 
w(A) — ^/XJM. Hence, the density of modes is 



piu:) = J d\p{X)SiL, - Lo{X)) ~ — , for F « J. 
Using the mean square displacement {x'^)u: ~ 1/Muj ^ F/cj, 



the spectral function results as, A{lo) « 



for w < F. Thus, these qualitative arguments are seen to re- 
produce correctly the Ohmic spectrum, its frequency range, 
and the F-independent coefficient of the replica solution (8). 
The latter is non-trivial given that both the mode density and 
the kinetic energy of the soft modes do depend on F. This non- 
trivial check makes us confident that the physical picture of a 
set of gapless, underdamped collective harmonic oscillators is 
the correct interpretation of the low energy spin excitations in 
this quantum glass. It is interesting to note that an analogous 
reasoning for spin glasses with metallic background leads to 
a similar picture, however with overdamped oscillators, and a 
spectral function growing as A{ijj) ^ l"^!^^^, again in agree- 
ment with replica theory [40, 52]. 

The appeal and potential of the TAP approach lies not only 
in the physical picture it provides for the collective excita- 
tions, but also in the fact that it lends itself to generalizations 
to finite dimensional glasses with long, but not infinite range 
interactions [56]. While in such more realistic models criti- 
cality might not exist at all length scales, collective modes are 
expected to persist down to very small energy scales, if the in- 
teraction range is large. Those may play an important role in 
activated transport of localized charge glasses [56], or induce 
non-Fermi liquid corrections in metallic glasses [57, 58]. It 
would be interesting to understand the spatial nature of these 
excitations, and to study whether and how the very low fre- 
quency spectrum eventually becomes dominated by more lo- 
calized, droplet excitations, as they occur in quantum Grif- 
fith phases. In the context of electronic glasses, many fur- 
ther interesting questions arise, in particular as to the interplay 
of glassy ordering with quantum phenomena such as (disor- 
dered) superfluidity and Bose glasses [5V], as well as Ander- 
son localization [ ]. 
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